December, 1996 
revised April, 1997 



IASSNS-HEP-96/130 
ITP_UH_28/96 
jr^jc/9612062| 



Scalar Deformations of Schwarzschild Holes 
and Their Stability 



00 
On 



Helge Dennhardt 2 and Olaf Lechtenfeld 1} 



1 School of Natural Sciences 
Institute for Advanced Study 
Olden Lane, Princeton, NJ 08540, U.S.A. 

00 

(N i and 

■ 2 Institut fur Theoretische Physik 

Universitdt Hannover 
^ \ Appelstrafie 2, 30167 Hannover, Germany 

CNJ ■ http://www.itp.uni-hannover.de/~lechtenf/ 

On 

cr 

5h . Abstract 

OX). 

We construct two solutions of the minimally coupled Einstein-scalar field equa- 
tions, representing regular deformations of Schwarzschild black holes by a self- 
interacting, static, scalar field. One solution features an exponentially decaying scalar 
field and a triple-well interaction potential; the other one is completely analytic and 
sprouts Coulomb-like scalar hair. Both evade the no-hair theorem by having par- 
tially negative potential, in conflict with the dominant energy condition. The linear 
perturbation theory around such backgrounds is developed in general, and yields 
stability criteria in terms of effective potentials for an analog Schrodinger problem. 
We can test for more than half of the perturbation modes, and our solutions prove 
to be stable against those. 
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1 Introduction 



It has been known for a long time that black holes are notoriously difficult to deform. This 
fact is captured by so-called no-hair theorems, which classify stationary, asymptotically flat, 
regular, black-hole solutions by a few conserved charges such as mass, angular momentum, 
electric and magnetic charges. Like any other no-go conjectures, no-hair theorems rest on certain 
assumptions, and it is worthwile to investigate those for necessity. 

The simplest case is that of "scalar hair" for Schwarzschild holes, where a gravitationally coupled 
static scalar field <1> is subject to some self-interaction potential V($). For a recent review, 
sec Q . The scalar no-hair theorem [^, [| rests on the dominant energy condition for the energy- 
momentum tensor T^, which stipulates that the energy current f = T^u 3 of the matter field 
should never be space-like. For a minimally coupled scalar field in a spacetime of signature 
( — h++), this translates to V(3>) > ||. Therefore, scalar deformations of Schwarzschild holes, 
if at all possible, need regions with negative interaction potential.^] 

The subject of this paper is to present two classes of static, isotropic, and regular solutions to the 
fully gravitating Einstein-scalar system, which circumvent the scalar no-hair theorem by virtue 
of having partially negative interaction potential. One of these solutions, with exponentially 
decaying scalar field, was given by one of the authors in an earlier paper Q. The second solution, 
featuring a Coulomb-like scalar field, is new and (in contrast to the first) completely analytical. 

After restating the starting point equations in Section 2, we discuss both classes of solutions 
in some detail throughout Section 3. In the second half of the paper, Section 4 provides a 
linear stability analysis of Schwarzschild hole deformations by general scalar fields. Following the 
techniques of Regge and Wheeler || , Vishveshwara Q , and Zerilli @ , we cast the perturbation 
equations into Schr6dinger-like form and read off stability criteria in terms of effective potentials. 
Section 5 finally specializes and applies the results to our two solutions, to decide upon their 
stability, before we conclude. 

2 The Equations 

The gravitational system we study in this paper is given by the Einstein-Hilbert-Klein-Gordon 
action 



for a self-interacting scalar field minimally coupled to standard gravity. The action is extremized 



We specialize to spherical coordinates x 1 = (t, r, (f). In the isotropic and static case all field 
degrees of freedom are functions of the radial coordinate r only, and the metric can be reduced 

1 The scale is set by the absence of a cosmological constant, V(&(oo)) — 0. 





by 



= V l di§-d$,V{<S>) , 



(2.2) 
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to two functions by residual coordinate transformations. The field configuration is then given 

$ = $(r) and ds 2 = -G(r) dt 2 + G{r)~ l dr 2 + R(r) 2 dtt 2 . (2.3) 

Instead of the usual Schwarzschild or isotropic coordinates, we have chosen the different but very 
useful gauge gttQrr = — 1. The equations ( |2.2[ ) reduce to three independent ordinary second-order 
differential equations, 

R" = 1& 2 R , 
4 

G» = ic*" + 2^1 , 

= \C¥ 2 + 2^ + 2^1fl , (2.4) 

where the prime denotes a derivative with respect to r. The equation of motion for $ follows 
from these. 



Before we construct solutions of ( p.4[) , it is necessary to put boundary conditions for small and 
large values of r. Near spatial infinity we take the metric to approach the Schwarzschild solution, 
i.e. 



G ^> 1-— and R ^> r , (2.5) 

r 



which, using ( |2.4D , implies that 

$ 1^2% $ Q + 0(r4) with v(<Z> ) = . (2.6) 

This precludes a cosmological constant. Linear stability of the asymptotic Schwarzschild metric 
yields the condition 

dlV($ ) >0 . (2.7) 

For small r, we like to encounter a black hole, signified by the presence of an event horizon at 
some value r = h where 

G{h) = and G(r > h) > . (2.8) 
The black-hole singularity at r = s is given by the (right-most) pole of G, 

G(r) ^> ±oo . (2.9) 

It proves useful to formally integrate the second equation of ( |2.4| ). With Q2.5| ), one arrives at 

oo 

, 2 f ~ 2f - 6M 



G(r) = R(r) 2 Jdr . (2.10) 

r 

The integration constant M is chosen to agree with the black- hole mass in ( |2.5| ) .0 Given a partial 
solution (&,R) of (|2.4|) , it remains to evaluate the integral above. 



2 -R(r) is not to be confused with the Ricci scalar. 

3 If $ decays as r -1 , then R ~ r + constant, and the two definitions of M differ by a constant. 
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3 Two Solutions 



In this section we will present two families of solutions to (2.4), one with Coulombic and one with 
exponential scalar field. The Coulombic solution is completely analytic and, to our knowledge, 
new to the literature. The exponential solution is partly analytic and was given earlier by one of 
the authors P|. If the potential is set to zero, the equations fl2.4| ) are soluble and yield (for fixed 
horizon) a known one-parameter family of solutions [||, |j, with $(r) ~ ln(l — h/r). Only for one 
value of the parameter one avoids the logarithmic singularity, at the price of setting = and 



recovering the Schwarzschild metric. For a general potential V($), the equations (2.4) cannot 
be solved analytically. We therefore approach the problem in reverse: 

• We make an ansatz for <£, choosing a class of functions 



We determine R from the first of ( |2.4| ). One integration constant is gauged away, the second 
one gets fixed by asymptotic flatness. 

Inserting R into ( |2.10| ), we obtain G. The mass M appears as an additional parameter. 
The right-most zero of G gives the event horion. 



• The third of (gj) yields U(r) = V($(r)). To find V($), we invert $(r) to r($). 

• The functions R(r), G(r), and V($) are plotted. 

3.1 Coulombic Scalar Field 

We assume that the scalar field has power-like behavior, 

= q-r a , (3.1) 
with the restriction a < — |, due to (|2.6| ). The resulting equation for R (2.4), 



4 

can be solved in terms of Bessel functions, 
R(r) = ^Z±_ (i— ■ r c 

2a \ 2 



R" = - q ^a 2 r 2a - 2 R , (3.2) 



2a \ 2 / 2a \ 2 



(3.3) 



The choice of sign originates in the $ — > — $ symmetry of the equations and will be hidden in 
the new coupling Q := ±q. 

In order to be able to proceed analytically, it is necessary to specialize further. Bessel functions 
become elementary for ^- = n + i, with n G Z, but the above restriction leaves only a = — 1, 
i.e. Coulombic fall-off for <£. Hence, we confine ourselves to the case 

$ = 2 . (3.4) 
r 
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Then, the solution ( |3.3D turns into 

r,/ \ / -2- -.fi\ -3- 

R{r) = C3 • r I e2r + e 2'- I + C4 • re 2r ; 

with new integration constants C3 and C4. We rescale time and radial coordinates by 

r = f/|c 4 | and t = i ■ | C4 1 
and absorb its effect on <£, G, R, and V by redefining 

g = q/\ci\ , M = M/\c 4 \ , c 3 = c 3 -|c 4 | 
so that C4 = ±1. Dropping the bars, R of (|3.5| ) (with C4 = ±1) tends to 

R > C3 • 2r ±r — * r , 

which allows only C3 = or C3 = 1, depending on the sign. In both cases, the result is 

R(r) = r ■ e 2r 

10.0 i — i — . 1 . 1 ■ 1 ■ ^7 . 1 
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r 

Figure 1: The function R(r) for three values of Q 

We must distinguish three cases: 

• Q > 0: R is a good radial coordinate, since R(r) is positive and monotonic. 

• Q = 0: The Schwarzschild case = 0). 

• Q < 0: R > everywhere, but R(r) is not monotonic. 

4 We may absorb a sign in the exponent into Q, since Q is denned only up to a sign. 



Inserting (|3.9| ) into (2.10) and evaluating the integral, one arrives at 
G(r) = 



2Q 3 

The asymptotic behavior, 



Q + 3M 2 Q + 3M 3M\ Q Q + 3M , 



Q 2 



G 



1 



2M - \Q 



signifies that the black-hole mass is actually 

M = 

in terms of which the metric function reads 



Gf(r) 



3Q + 6M 2 3Q + 6M Q + 6M\ q 3Q + 6M 



4Q3 



r + 



2Q 2 



• r 



2Q J 



e<~ + 



4Q3 
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a 
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Figure 2: The function G(r) for M = 1 and three values of Q 



Again, a case distinction should be made: 

• Q > 0: the scalar field increases the size of the black hole. 

• Q = 0: the Schwarzschild solution. 

• — 2M < Q < 0: the scalar field decreases the size of the black hole. 



Q < —2M: a naked singularity, since G ► +oo without a horizon. 



Therefore, we restrict ourselves to 

Q > -2M 



(3.14) 



Finally, we employ the third of fl2.4|) to compute the potential, 



U(r) 



3(Q + 2M) 

Q 3 



Q 



3 H — 7T • sinh — — 3 — • cosh — 



Q 



and invert via r = q/& to obtain 



V(<&) 



(3 + $ 2 ) • sinh |^| - 3|$| • cosh<& 



with 



A 



3(Q + 2M) 



> 



(3.15) 



(3.16) 



(3.17) 



Since $ = g/r reaches only positive (or only negative) values of <£, (3.16) is valid only for one 
sign of <£. It is, however, natural to impose a $ — > — $ symmetry, and ( 3.16 ) respects it already. 
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Figure 3: The function V($) for A = 1 



The potential function V(3>) is not analytic at the origin, but varies as 



V(») 



15 1 1 



(3.18) 



so that d^V(0) = and the vacuum stability is marginal. It is apparent that the interaction 
potential is everywhere negative. Therefore, the no-hair theorem |3) simply does not apply to 
our solution since the latter does not obey the dominant energy condition. 
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3.2 Exponential Scalar Field 

We assume that the scalar field drops exponentially, 

= q-e- mr , (3.19) 

with m > for asymptotic flatness. This ansatz and the solution presented in the following were 
the subject of Q. 

Again, the resulting equation for R, 

R"(r) = ^q 2 m 2 e- 2mr R(r) , (3.20) 
can be solved in terms of Bessel functions, namely 

R(r) = Z U^- ■ e"" 

= Cl • Jo ■ e~ m A +C2 -K (^- e~ mr ^ , (3.21) 

with two integration constants c\ and c<i. Here too, we rescale time and radial coordinates by 

r = f/m and t = i ■ m (3.22) 
and absorb its effect on <£, G, R, and V by redefining 

M = M/m (3.23) 
so that m = 1. With a second coordinate transformation (7 is Euler's constant), 

r = r + In 



+ 7"- , (3.24) 

C2 



we get rid of c\ , but must redefine 



M = M + - ■ { In 

3 



+7 — 

C2 



4 

q = sgn(g) -4e Cl/c2 - 7 . (3.25) 
Dropping the hats, R of ( p.21|) (with m = 1 and c\ = 0) tends to 

R c 2 -r -U r , (3.26) 

so that we have to take C2 = 1. Inserting the infinite series representation of Ko, the radial 
function may be expressed as 

R(r) = r + £( r + £i)-^ ? (fe-f . (3.27) 

k=i e=o y ' 

The sign ambiguity of ( |3.21| ) has disappeared since only q 2 enters. 



1.00 



0.75 - q=3 




r 

Figure 4: The function R(r) for four values of q 



Three cases are to be distinguished: 



\q\ > 4exp(— 7): R > everywhere, but R(r) is not monotonic. 

< \q\ < 4exp(— 7): R is a good radial coordinate, since R(r) is positive and monotonic. 
\q\ = 0: The Schwarzschild case (3> = 0). 



Inserting (|3.27|) into ( glo| ), one discovers that the integral cannot be done analytically. It is 
sensible, however, to expand the integrand in powers of q/4 ■ exp(— r) and integrate it term by 
term. The result is the series 



00 00 
G(r) = Y,[ a k(r,M)+Y,he(r,M)-q 2l Ei(2£r) 



k=0 



2k -2kr 

1 & 



(3.28) 



where the coefficients a& and bki are linear in M and polynomial in r and -. Ei(x) = f^°dy e~ y /y 
denotes the exponential-integral function. Asymptotically, ( p. 28 ) yields the Schwarzschild metric 
with mass M. The following graphs were obtained by keeping terms up to k = I = 6, i.e. 
0(exp(-12r)). 

Despite the fact that R depends greatly on q (see Figure |I|) , the metric function G hardly varies 
with q. Magnifying Figure |5| brings out a small shift of the horizon. Note that, contrary to the 
Coulombic case, the singularity is shielded for all values of q. 

We finally come to the potential. Carrying over the series expansion to the third of ( p.4|) , one 
gets an expression of the same form as (3.28), 



00 00 
U(r) = Y,[^(r,M)+Y J d k e(r,M)-q 2e Ei(2£r) 



. q 2k e -2kr 



(3.29) 



fc=i 
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Figure 5: The function G(r) for M = 1 and four values of q 
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Figure 6: Detail of Figure |5| near the horizon 



where and d^e share the properties of a& and 6^. Inverting by r = — In \ $/q\ we find 



oo oo 

V@) = 5][4(l^l>M)+^4,(|$|,M)-^Ei(-2€ln|$/g|) 



2k 



(3.30) 



fc=i 



where c' fc (|<l>(r)|, M) = Cfc(r, M) and d' w (|$(r)|, M) = du(r,M). Again, it is natural to extend 
the solution for only positive (or only negative) values of <1> to the whole real axis by taking 
V{—$) = V{$), already made explicit in (|3.30|). 
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q=3 
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Phi 

Figure 7: V(<&) for M = 1 and three values of q 
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Figure 8: V(3>) for q = 1 and three values of M 



Taking a closer look at the origin reveals the local minimum which also follows from 

V($) > ci • $ = -$ - - r^T-i 

v ; 2 In $/g 



(3.31) 



and, again, is nonanalytic but smooth (even C°°). So, contrary to the Coulombic solution, 
vacuum stability is manifest. 

The no-hair theorem [||, [3| gets circumvented by having regions of negative potential, again vio- 
lating the dominant energy condition. In 0] , the energy-momentum tensor of this configuration 
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0.0005 




was displayed. Near the horizon, the energy density p = V + \G& drops below zero, which 
violates the weak energy condition as well. 

4 Perturbations 

To examine the stability of the static solutions introduced in Section 3, we first decompose 
possible perturbations into tensorial modes and simplify them with a gauge transformation and 
a rotation. We then use this general ansatz for the perturbations to linearize the equations 
of motion around the static solution and put them into a Schrddinger-like form. In this form, 
instability is equivalent to the existence of negative eigenvalues of the Schrodinger problem, a 
fact which leads to a general stability condition for static solutions. 

This procedure for dealing with stability questions in general relativity was first described in 
1957 by Regge and Wheeler ||, who proved the stability of black holes against odd parity 
perturbations. The proof for even-parity perturbations was achieved in 1969 by Vishveshwara 
|]], using Kruskal coordinates. The work of Zerilli |7| completed these examinations. 

In this paper, we are dealing with a more general system containing an additional scalar field, 
which implies that it is always possible to get the equations of Regge, Wheeler, Vishveshwara, 
and Zerilli by carrying out the zero scalar field limit ($ -> 0, F(*) -> 0, G -> 1 - R -> r). 
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4.1 Canonical Ansatz 



To get a general ansatz for the perturbations of our system, we first separate the perturbations 
from the background 



with 5gij and <5<J> being infinitesimal small. 

Described in this way, the perturbations are independent fields, living on a static, spherically 
symmetric background, which implies constancy of energy, angular momentum, and parity. Be- 
cause of this, modes of different energy, angular momentum, or parity do not mix when time 
evolves, and each mode can be examined separately for stability. 

The decomposition into modes with fixed energy is achieved by introducing a factor exp(iwi) in 
each mode. Likewise, the decomposition into modes with fixed angular momentum yields factors 

Yff{p,<p)- This leaves 

a pure radial part for the remaining degrees of freedom, which then is 
split into two parts with opposite parity. 

The functions Yff already have definite parity (— 1) L , so the parity separation is not necessary 
for the scalar perturbations 5&. In contrast, the metric perturbations 5gij with their tenso- 
rial character bear directional information which does not remain unchanged under a parity 
transformation and hence must be split carefully into two parts. This is shown i.e. in pfl. 

The resulting formulas can then be simplified enormously by a gauge fixing, which eliminates 
four radial degrees of freedom. Furthermore, a rotation is performed for each mode, which sets 
M to zero and thus transforms the functions Yff into Legendre polynomials Pl(cos^) so that 
all dependence on tp disappears. Regge and Wheeler called the resulting perturbation modes 
canonical^]. They read 

for even parity: 




9ij + Sgij 





V o 






hj(r) sin 2 t?/ 







for odd parity: 




/ h 8 (r)\ 

h 9 (r) 



\h 8 (r) h 9 (r) / 




<5$ 







(4.3) 



where "even" and "odd" denote parity (— 1) L and (— 1) L+1 , respectively. 



5 To get the notation of 0|, replace hi — » Ho ■ G, hi — > Ho, /13 — > H2/G, hr — > K ■ R 2 , h% — * ho, hg — » h 
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4.2 Odd-Parity Sector 



To examine the odd-parity perturbations, we first insert Q4.1|) with the canonical odd-parity 



ansatz (4.3) into the equations of motion ( p.2|) and retain only terms linear in the perturbations 
hg and hg. 

Along this procedure, one encounters second and third order derivatives of the Legrende poly- 
nomials, which can be reduced in order by their defining property 



d|P L (costf) 



COS 

sini? 



dtfP L (costf) - L(L + 1) • P L (cos i?) 



(4.4) 



so that only derivatives of order and 1 are left. After tedious calculations one obtains two 
independent differential equations, 



+ i)-i i 



R 2 

sin^^P L (cos^)e i£J * 



R! 



1 iu) 



2G)- h9 + lUJ GR- h8 -2G' h ' 8 



(4.5) 



— ■ h 8 - G' ■ hg - G ■ tig 



^L(L + 1) sin P L (cos 0) + cos d,#P L (cos 0) ) • e iujt 



(4.6) 



These are the elements (1,3) and (2,3) of the perturbation equations for the gravitational field. 
Equation Q4.5|) was simplified by using the static equations of motion (2.4) to eliminate R" and 
V. 

The element (0,3) of the graviational equations results in a dependent differential equation, all 
other elements are given by symmetry or yield the identity. The equation for the scalar field also 
yields nothing, since the scalar perturbations do not take part in the odd-parity oscillations at 
all. The system is hence fully described by the two differential equations (|4.5| ) and fl4.6|) . 

At this point, it is necessary to separate out two degenerate cases for low values of the angular 
momentum L. They are 

• L = 0: In this case, Po(cos0) = 1, both differential equations (|4.5| ) and ( |4.6| ) degenerate 
to the identity, and no odd-parity perturbation really exists, since they vanish totally in 



our ansatz (4.3) 



L = 1: In this case, Pi (cos 0) = cost? and equation ( |4.6j ) degenerates. An additional gauge 
transformation, 



x 



x l + 5'' 



*uR 2 



ti ■ e 



(4.7) 



and a redefinition of hg make hs vanish totally from our perturbation ansatz. The remaining 
differential equation ( [4.5[) then states hg to be zero as well. Hence, the perturbation has 
been gauged away completely. 



13 



This discussion shows that the lowest angular momentum of odd-parity perturbations is L = 2, 
hence the angular parts of the differential equations never vanish. Since these equations have 
to hold for all angles and all times, the radial parts (in square brackets) can be set to zero 
separately, leaving a pure radial differential system. 

Solving the radial part of (p~6|) for h%, one gets 



Inserting this into Q4.5Q yields 



h 8 = — ■ (G' ■ hg + Gh'o) . (4.8) 



(2^ - G' 2 - GG" - £ (L(L + 1) - 2)) • hg 

+ (2^^ - 3GG'^J -h' g -G 2 - h'l = u? ■ hg . (4.9) 



As one can see, equation (|4.9| ) determines the entire dynamics, while (4.8) only expresses h% 
in terms of the dynamical variable hg. For vanishing interaction potential, it reduces to the 
Regge- Wheeler equation Q. The remaining problem now consists of solving the second-order 
differential equation ( |4,9| ), which in fact is not possible without specifying the static solution 
functions G and R. 

For stability questions, however, it is unnecessary to solve this differential equation explicitely. 
We only need information about the spectrum of frequencies u. Because of the time evolution 
exp(iujt) of the perturbation modes, it is sufficient to determine whether uj has an imaginary part 
or not. To decide this question, we will put our differential equation ( [4.9D into a Schrddinger-like 
form. To this end, we introduce a new dynamical variable ^ via 

Hr) ■= |^-Mr) (4.10) 

and also the tortoise coordinate x by 

dx := — --dr . (4.11) 

Since G > everywhere outside the horizon, the new coordinate x is a regular, monotonic 
function of r there. Thus, it may serve as an alternative radial coordinate for the exterior of 
the black hole. Because of asymptotic flatness, G comes close to 1 for large distances, and thus 
x approaches infinity there. Near the horizon, G becomes zero, and hence x tends to negative 
infinity. So x parametrizes the whole exterior of the black hole. Let us also look at the new 
dynamical variable ^f. It is regular everywhere outside of the horizon, as G and R are regular 
there and R > 0. For x — > +00, \I/ approaches zero faster than 1/x, but for x — ► —00, * decreases 
even exponentially. Thus, \P turns out to be a Fourier-transformable function. 

With these new features, equation Q4.9p now reads 



-fi%V(x) + V e ff(x)V(x) = u) 2 ^{x) , with (4.12) 

v '"^ = §B iL(L+i) - 2)+Rix) { a 'W)) • <4J3) 
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which in fact is in Schrodinger-like form. The effective potential belongs to a special class of 
potentials; it is regular everywhere and tends to zero at both infinities. For this kind of potential 
in one-dimensional Schrodinger problems, it is well-known that there may exist a bound state 
only if the potential is negative somewhere (see e.g Q). A bound state corresponds to a negative 
eigenvalue w 2 , i.e. an unstable mode must exist. So a sufficient condition for stability is 

i 

V ef f(x) > \/x G [—00, +00] . (4.14) 

A further look at ( t4.12| ) reveals an angular momentum barrier j^(L(L+l) — 2), which grows 
monotonically with L, so that it is sufficient to perform the test on the mode with smallest value 
of L. The smallest non-trivial L equals 2, and thus the final stability condition, written in the 
old coordinate r, reads 

! 

V e ff(r) > Vr€[/i,+oo] , with (4.15) 

V < ^ i G ^ 1 M r ) 2 R'( r f G{r)G'(r)R'{r) G(rfR"(r) 
Veff{T) ~ 4 i?(r) 2 + R(r) 2 R(r) R(r) ' 

4.3 Even-Parity Sector 

Finally, we examine the even-parity perturbations. Again we first insert fl4.1|) , this time with 
the canonical even-parity ansatz ( |4.2| ), into the equations of motion ( |2,2D and keep only terms 
linear in h\,h2,h^,hi, and (p. Again one encounters derivatives of Legendre polynomials and 
reduces them via their defining properties. The resulting differential system is very complicated 
and difficult to handle. It seems impossible to achieve a general Schrodinger-like form in the 
canonical gauge, as we just did in the odd-parity sector. We are searching for an alternative 
gauge and work on this subject is in progress. 

Nevertheless, it is possible to achieve a Schrodinger-like form of the differential equations in the 
special case L = 0, which represents monopole oscillations. For L = it is possible to gauge 
away two further degrees of freedom and we choose 

hi = and h 7 = . (4.17) 

This additional gauge-fixing results in an enormous simplification of the differential equations. 
The (0,1) component of the gravitational perturbation equations is now easily solved for /13, 

1 R& 

Given this result, it is possible to solve the combination (2,2)+(2,3)/ sin 2 $ of the components 
of the perturbation equations for h,2, 

i , RGQ" R 2 G& 3 , i GR& ,, 



We have simplified the equation by using the static equations (|2.4f ) to eliminate the function 
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Taken together, this allows us to decouple the scalar perturbation equation for (f> ( |4,2| ) from the 
metric perturbations, 

/ 1&*R 2 G 3__, 2 &RG&' 1& 2 RG' u 2 
+ (2^ + G'j ■ 4>' + G ■ 0" . (4.20) 



R 

What is more, all other components of the perturbation equations turn out to be dependent 
or degenerate to the identity and, hence, give no further information. It is obvious that the 
second-order differential equation ( 4.20| ) determines the whole dynamics of the system, while 



(gig ) and (gig ) only express /12 and I13 in terms of the dynamical variable 



As in the odd-parity sector, it is possible to put the remaining equation (4.20) into a Schrodinger- 



like form by introducing a new dynamical variable $ and the tortoise coordinate x, 

(r) := R(r) ■ <f>(r) , dx := —— ■ dr . (4.21) 

G{r) 

The coordinate x again parametrizes the whole exterior of the black hole, and the new variable 
is again regular everywhere. Because of asymptotic flatness, R —* r for large distances, and 4> 
like 1/r or faster. This means that \& approaches a constant at spatial infinity. The same behavior 
can be seen at the horizon, where x tends to negative infinity. approaches another constant 
here, so that ^ turns out to be a finite function and thus Fourier-transformable. Expressed in 



the new variables, equation (4.20) now reads 



-d 2 -^ + V e ff^ = uj 2 -^ , with (4.22) 

V,U ^ ^ - l** - ^ - \^ + ^ + , (4,3, 

and primes denote derivatives with respect to r. The effective potential V e // shows a slightly 
different behavior at the two boundaries, namely 

lim V eff (x) = lim V eff (x) = a|V(* ) • (4-24) 

There are two possibilities: 

• d\V{<&o) < 0: It is always possible to find a solution for ^ with negative eigenvalue co 2 , 
thus the system is unstable. 

• c^y^o) > 0: A positive semidefinite effective potential is sufficient but not necessary 
to exclude bound states. If it is negative somewhere, a further investigation, perhaps of 
numerical nature, must be performed. 



5 Stability 

The results of the previous two subsections suffice to examine the stability properties of the static 
solutions constructed in Section 3. We shall apply the general criteria based on the effective 
potentials ( 4.16| ) and (4.23) to each solution in turn. 
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5.1 Stability of the Coulombic Scalar Field Solution 



We first look at the odd-parity sector of the static solution with power-like decaying hair (Sec- 

The effective 



10 



tion 3.1). The effective potential ( 4.16 ) for this sector is shown in Figure 
potential vanishes exactly at the horizon (whose location moves with Q) with a positive slope, 
and falls to zero exponentially at large distance. It is positive everywhere exterior to the black 
hole, signalling stability. 

For the even-parity monopole sector (L = 0), the second of the above two cases is relevant, 



since d$V(0) = for the Coulombic solution (|3.16 ). Thus, the monopole sector is stable if the 



effective potential never turns negative. The effective potential in this sector is shown in Figures 



11 and 12. By inspection, the potential becomes negative for too large or too small values for 
the parameter Q, so that the system may become unstable there. The upper stability bound for 
Q can be found algebraically. Demanding positivity of the leading term at large distances, one 
finds 



Q 



< -M 
2 



(5.1) 



The lower bound seems only accessible numerically. For M = 1, one gets roughly Q ~ —0.8. 

To summarize, the Coulombic scalar field solution (Section 3.1) is completely stable against 
odd-parity perturbations, and is stable against even-parity monopole perturbations at least for 
a certain range of Q. The stability in the L>0 even-parity sector could not be decided here. 
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Figure 10: V e ff of the odd-parity sector for M = 1 and several Q > —2M 



17 




18 



5.2 Stability of the Exponential Scalar Field Solution 



As for the other static solution (Section 3.2), we first look again at the odd-parity sector. The 
effective potential ( |4.16| ) for this sector is shown in Figure It is surprising to observe that 
the effective potential changes only minimally with q, even though the static solution function 
R varies dramatically within the same parameter domain (compare with Figure ||). A zoom 
(Figure |l4| ) reveals that the effective potential does vary with q, as its zero always sits at the 
(moving) horizon. The effective potential never turns negative, so the system is stable against 
odd-parity perturbations. 

Finally, we turn to the stability in the even-parity monopole sector. This time, d$V(0) = 1 > 
from ( 3.30 ), so again we deal with the second case above. This means that positivity of the 
effective potential is a sufficient but not necessary criterion for stability. The effective potential 
for the monopole sector is shown in Figure 

H As one can see there again, the effective potential 
depends only weakly on the parameter q, while the static solution, especially R, changes dra- 
matically. A closer look at the horizon (Figure |i~6| ) shows the same feature as in the odd-parity 
sector. With the effective potential being positive in the complete exterior of the black hole, the 
system is stable against monopole perturbations. 

To summarize, the exponentially decaying scalar field solution (Section 3.2) is stable both against 
odd-parity and monopole perturbations. The answer to the stability question for higher angular 
momentum in the even-parity sector still eludes us. 
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Figure 13: V e ff of the odd-parity sector with M = 1 and several q 
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Figure 14: Zoom of the odd-parity V e ff around the horizon with M = 1 
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Figure 15: V e ff of the monopole sector with M = 1 and several 
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Figure 16: Zoom of the monopole V e ff around the horizon with M 
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6 Conclusions 



We have demonstrated that it is possible to obtain asymptotically flat regular black-hole solu- 
tions of the static and isotropic Einstein-scalar field equations, if one gives up the dominant 
energy condition. These deformed Schwarzschild holes avoid the no-hair theorem by having par- 
tially negative scalar self-interaction potential. In some cases, the metric and scalar configuration 
can be computed analytically (we gave one example), and a generic feature seems to be a mild 
non-analyticity of the potential at zero field. 

The indefinite scalar potential does not imply instability however, at least not against general 
odd-parity or monopole even-parity perturbations. For these, we could derive sufficient stability 
conditions, which we then tested on our examples, with positive results (in some parameter 
range). Of course, this does not yet preclude potential instability against even-parity perturba- 
tions with higher angular momentum. Indeed, a numerical investigation seems to point in that 
direction |lC[| . 

There are a number of obvious extensions and open questions, in particular 

• Are there more analytic solutions? 

• Can one add a cosmological constant? 

• What effect does the scalar field have inside the horizon? How do our solutions look in 
Kruskal coordinates? 

• With an additional gauge field, can one still find (partially) analytic solutions? 

• May one generalize our solutions to the axisymmetric situation, using the Ernst potential? 

• Is there a minimal amount of violating the dominant/weak energy condition for hair? 

Some of these issues are under investigation, and we hope to report on them in the near future. 
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